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Abstract. - In an elegant paper, D. Nelson suggested a method to produce tetravalent 
colloids based on a tetrahedral configuration created on the surface of a spherical particle. It 
emerges from a two-dimensional nematic liquid crystal placed on a sphere due to the presence 
of four 1/2 disclinations, i.e., topological defects in the orientational order [1]. In this paper we 
show that such a tetrahedral configuration also occurs in the wetting layers which form around 
spheres dispersed in a liquid crystal above the nematic-isotropic phase transition. Nematic 
wetting therefore offers an alternative route towards tetravalent colloids. 



Recently, D. Nelson proposed a method by which tetravalent colloids could be produced [1]. 
These micron-sized particles would have four chemical linkers or DNA strands symmetrically 
attached to their surfaces. Similar to, e.g., carbon, silicon and germanium atoms with their sp 3 
hybridized chemical bonds, the tetravalent colloids could then arrange into, e.g., a colloidal 
crystal with a diamond lattice structure which is known to possess a large photonic band 
gap [2]. 

To create the attachment sites for the chemical linkers, Nelson considered the two-dimen- 
sional nematic liquid crystal phase on the particle surface realized by a monolayer of elongated 
constituents such as gemini lipids [3] , ABA triblock copolymers [4] or nanorods [5] . One could 
first think that the rodlike particles on a sphere create Mermin's boojums [6], i.e., point defects 
with topological charge 1 situated at the north and south pole of a sphere as illustrated in Fig. 

The lines indicate the director field. A vector order parameter would indeed create such 
a configuration. The nematic order, however, is described by an axis in space which allows the 
boojums to split up into a pair of disclinations with charge 1/2. Lubcnsky and Prost showed 
that the ground state of a 2D nematic texture on a sphere consists of four 1/2 disclinations 
situated at the vertices of a tetrahedron [7] . Fig. ^)) tries to illustrate such a configuration 
with tetrahedral symmetry in the angular space of spherical coordinates, i.e., polar angle 9 
versus azimuthal angle 4>. For the moment, let the rectangles symbolize the nematic director. 
To visualize the boojum configuration in this representation, all the local directors would 
point along the vertical (lines of longitude). Squares situated on the 9 = 0,7r line (north and 
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Fig. 1 - Different configurations of a nematic liquid crystal on a sphere. The pictures apply to both 
a pure two-dimensional nematic and to biaxial nematic wetting layers, a) Boojum configuration with 
two +1 disclinations at the north and south pole. The lines indicate the director field, b) Tetrahedral 
configuration visualized in the angular space of spherical coordinates. The orientational order in the 
sphere is indicated by rectangles. Four + 1/2 disclinations (see dots) are located on the vertices of a 
tetrahedron. As a reference, the square in the upper left corner symbolizes a fully isotropic order of 
the molecular axis. 

south pole) would indicate the isotropic order in the core of the boojums. To arrive at the 
tetrahedral configuration, the boojums split up into a total of four 1/2 disclinations which 
move away from the north and south pole as indicated by the dots in Fig.^). If one is able to 
attach chemical linkers selectively to the core of the disclinations, then a tetravalent colloidal 
chemistry can be realized [1]. 

Does the tetrahedral configuration survive when the colloid is suspended in the nematic 
phase of a three-dimensional liquid crystal? The 1/2 disclinations would extend as lines into 
the solvent. For topological reasons, they have to end either on neighboring particles or at 
the boundaries of the system. However, the resulting network of defect lines carries a lot of 
free energy. It can be greatly reduced if the particles realize the boojum configuration. The 
corresponding +1 defect lines in the bulk nematic phase are not stable and "escape into the 
third dimension" [8] . This is indeed observed in nematic emulsions [9] . 

We have recently looked at nematic wetting transitions and capillary condensation in 
connection with colloids suspended in a liquid crystal solvent above the nematic-isotropic 
phase transition [10]. So the interesting question arises if the tetrahedral configuration occurs 
in nematic wetting layers which surround colloidal particles. This would then suggest an 
alternative route towards tetravalent colloids. In answering the question, one first needs to 
establish orientational order close to the colloid with a preferred axis parallel to its surface. 
Since the orientation of the preferred axis should be free to rotate within the surface we cannot 
just fix it by an appropriate surface potential. All we can do is to favor planar ordering of 
the liquid crystal molecules at the surface with an isotropic distribution of their axes in the 
surface. This uniaxial oblate order is, however, energetically disfavored in the bulk since close 
to the phase transition the molecules already have a tendency to align parallel to each other, 
i. e, to assume uniaxial prolate order. As a compromise, biaxial order close to the bounding 
surfaces occurs and introduces a preferred axis as required. The phase transition from uniaxial 
oblate to biaxial surface ordering is well studied in literature for planar geometry [11-13]. In 
the following, we readdress this problem to set the stage for the spherical geometry. For the 
latter we then show that close to the nematic-isotropic phase transition, the biaxial ordering 
indeed leads to the tetrahedral configuration as in the two-dimensional case. In addition, we 
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find that the transition from oblate ordering to the tetrahedral texture occurs via the boojum 
configuration which possesses a narrow stability region only. 

To quantify the surface-induced orientational order, we use a traceless and symmetric 
second-rank tensor Qij, also called alignment tensor. We perform our analysis with the help 
of the Landau-Ginzburg-de Gennes free energy density in the bulk [14], f(Q tj ) = |a (T - 
T*)QijQij - bQijQj k Q k i + c(QijQij) 2 + |-ki(Qij,jfe) 2 + \L 2 (Qi ] .j) 2 where summation over 
repeated indices is implied and the symbol , k means spatial derivative with respect to Xk ■ The 
first three terms describe the nematic-isotropic phase transition; do and c are positive constants 
and T* denotes the supercooling temperature of the isotropic phase. The fourth and fifth term 
penalize any non-uniform orientational order. Since we are interested in the basic features of 
our system, we always choose L\ = L 2 - The anchoring of the molecules to a bounding surface 
is quantified by the Nobili-Durand free energy density, fs(Qij) = ¥f(Qij — Qij'') 2 ! where W is 
the anchoring strength and the preferred order parameter at the surface [15]. The number 

of parameters is reduced by using a rescaled order parameter /iy = Qij/ s [s = 6/(v6c)] and 
temperature r = 12ceto(T — T*)/b 2 . Furthermore, all lengths are given in terms of the particle 
radius a (r — r/a) and the unit of the free energy is A/a 3 = b 4 a 3 / (36c 3 ). We also introduce 
the nematic coherence length at the nematic-isotropic phase transition, £ r = (12cLi/6 2 ) 1//2 , 
and define the reduced anchoring strength as 7 = 6cW/(b 2 £ r ) so that the reduced free energy 
becomes 

+ik J *, . M 

Note that in the planar case, we choose a = £ r as the unit length. Typical values for the 
nematic compound 5CB are Af = 8 • 10 5 erg/cm 3 , £ r = 10 nm, and the temperature interval 
At = 1 corresponds to 1.12 K. 

A uniaxial order parameter /iy = ^S(niiij — Oy/3), where n is the nematic director and 
Sij the Kronecker symbol, minimizes the first three terms on the right-hand side of Eq. (I). 
The bulk nematic-isotropic phase transition from S = to Sb = v6/4 occurs at r c = 1, and 

=9/8 is the superheating temperature of the nematic phase. The bulk exhibits prolate 
order since Sb > 0, i.e., the rodlike molecules align along a common axis. In the following, 
we assume that bounding surfaces favor oblate order. The preferred order parameter at the 
surface therefore assumes the form /z^ = | < S'o(i^£'j — Sij/3) with So < 0. In general, the 
order parameter can be biaxial. We visualize it by adding a term proportional to 5ij and by 
calculating its eigenvectors and eigenvalues whose directions and magnitudes then define the 
edges of a cuboid. We find that in the planar or spherical geometry one eigenvector always 
points along the surface normal or the radial direction, respectively. So in Fig. ^>) we only 
see one face of the cuboid. To test for the biaxiality of /x^-, we introduce the measure [16] 

R 1 a (^ij^jk^ki) 2 

which is zero for uniaxial tensors and assumes a maximal value of one. An averaged degree 
of biaxiality is defined as 

- = Jd?fHijtJ.ijB ^ 

It weights the local B with the degree of orientational order given by (J>ij(J>ij such that for 
fiij — > unphysically high local biaxialities do not contribute to B. 
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Stable orientational textures correspond to a minimum of the free energy functional ■ Its 
variation gives Euler-Lagrange equations for the five independent components of fiij (f). They 
are discretized and then solved numerically by the standard Newton-Gauss-Seidel relaxation 
method using appropriate starting configurations [17]. The planar geometry is just a one- 
dimensional problem with the bounding surface at z = 0. To deal with the infinite half space, 
we introduced a new coordinate p = exp(— z). Since close to the bulk phase transition the 
biaxial wetting layer exhibits complete wetting [12], we also used an adaptive grid to resolve 
the interface between isotropic and orientational order. Oblate and biaxial wetting layers were 
distinguished by the biaxiality parameter B at the surface. The wetting layers close to a sphere 
were determined in a spherical coordinate system using the local coordinate basis to define 
the alignment tensor ^Ar). We replaced the radial coordinate f by p = exp[— (f — l)a/£ r ] 
to account for the infinite space around the sphere. The Euler-Lagrange equations were 
formulated with the help of differential geometry. Phase transitions between different wetting 
layers were monitored by the averaged biaxiality B and by the behavior of the free energy 
which we calculated numerically by applying Simpson's integration rule to Eq. |T|. 

We add some comments about the possible orientational textures close to a sphere. It has 
already been noted that one eigenvector of the tensor order parameter always points along the 
radial direction. From the remaining eigenvectors, we take the one with the largest eigenvalue 
as a director. Close to the particle it can then form the splay dominated boojum of Fig. 
or the bend dominated version where the director always points along the lines of latitude 
of the sphere. Within our numerical accuracy, we observe that both configurations have the 
same free energy. This is due to the fact that in the Landau-De Gennes theory splay and 
bend deformations are energetically degenerate. The same applies to the bend dominated 
tetrahedral configuration of Fig. \\]p) and its splay dominated counterpart which follows from 
Fig. by rotating all rectangles by 90°. Note that Fig. QJ}) illustrates the orientational 
order a distance of 0.05a away from the particle surface. We found that around this distance 
the biaxiality of the orientational order is the largest. Topological defects in biaxial nematics 
are more complex than in the uniaxial case [18]. However, since one eigenvector always points 
along the radial direction, we do not have to worry about these complications and can basically 
work with the director introduced above to identify the defects. 

Before continuing the discussion of the spherical case, we summarize our results for the 
planar geometry in Fig. [3 We show a phase diagram in terms of temperature r versus reduced 
anchoring strength 7, the preferred surface order parameter So is fixed for each transition line. 
Above Sq ~ — 1 only the uniaxial oblate phase exists. For smaller So, a biaxial phase appears 
at temperatures close to the bulk phase transition (r = 1) and for sufficiently large anchoring 
strength. Decreasing So widens the existence region of the biaxial phase. The transition 
from the uniaxial oblate to the biaxial wetting layer is either of first (solid line) or second 
(dashed line) order. The transition lines meet in a tricritical point (dot). For large 7, the 
lines approach a constant temperature r > 1 (note the logarithmic scale for 7). Our results 
on a single transition line agree with earlier findings where a one-parameter surface potential 
with just a surface ordering field was used [12]. With our more realistic surface potential, we 
show in addition that the biaxial wetting layer only appears for sufficently small So- 

Based on the knowledge of the planar case, we now address the spherical geometry. For 
300 nm particles (a/£ r = 30), Fig. presents a three-dimensional phase diagram in terms 
of temperature r, reduced anchoring strength 7, and preferred order paramter Sq locating 
the different types of wetting layers. The phase diagram is the result of extensive numerical 
calculations which determined the three-dimensional order parameter field close to the particle. 
They were performed on a PC cluster and would take a few years on a typical state-of-the art 
PC. As in the planar case, the uniaxial oblate wetting layer is always stable for < So < — 1. 
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Fig. 2 - Phase diagram for the planar geometry in terms of temperature r and reduced anchoring 
strength 7; the parameter of the curves is the preferred surface order parameter So. First and second- 
order phase transitions between uniaxial oblate and biaxial wetting layers are indicated by solid and 
dashed lines, respectively. The dot means tricritical point. 



Decreasing Sq induces first a phase transition to a boojum configuration for sufficiently large 
7 and then to the tetrahedral wetting layer. Note that the stability region of the boojum is 
restricted to a very narrow So intervall and, therefore, hardly visible in Fig.[3Ji). The solid line 
located on the upper surface indicates a line of tricritical points. To the left of this line, the 
transition from oblate ordering to the boojum is of first order and to the right it is of second 
order. In contrast, the transition from the boojum to the tetrahedral configuration is always 
discontinuous. We never found that the boojum evolves continuously into the tetrahedral 
configuration via an asymmetric arrangement of the 1/2 disclinations. One could speculate 
about the narrow stability region of the boojum configuration. After all, as noted above, in 
the bulk nematic phase the corresponding +1 disclination lines are not stable and "escape 
into the third dimension" to reduce the free energy. This, however, cannot occur in our case 
since for the biaxial order parameter the equivalent +1 disclination line is actually stable [18]. 
So the system reduces the free energy by adopting the tetrahedral configuration. 

In Fig. Eb) we draw contour lines of the surfaces of the three-dimensional phase diagram 
for different values of So using the spline function of the graphical program gnuplot. The 
resulting two-dimensional phase diagrams resemble very much the one in the planar case (see 
Fig.[2J) with the exception that for decreasing temperature, we have a sequence of three phases 
as indicated for So = —2. As before, the solid and dashed lines mean first and second-order 
transitions, respectively, and the dots show the tricritical points. To account for size effects, 
we also studied particles with radius a = 100 nm (a/£ r = 10). The phase diagrams in Fig. 
0] qualtitatively look the same as for the larger particles. Comparing the two-dimensional 
diagrams of the planar case and the two particle sizes for a fixed So , one clearly recognizes 
that with increasing curvature the temperature range where the biaxial wetting layers exist 
becomes smaller. So curvature suppresses the biaxial phases. Surprisingly, the rough location 
of the transition lines with respect to the anchoring strength 7 does not seem to be affected 
by curvature. 

How reasonable are the rescaled parameters of the surface potential where the biaxial 
wetting layers are observed? Using Landau parameters from 5CB (b — 0.7 ■ 10 7 erg/cm 3 , 
c = 0.5 • 10 7 erg/cm 3 ) [19] and £ r = 10 nm, we arrive at 7 = W/l-^j. So the well achievable 
anchoring strength W — lerg/cm 3 gives the right order of magnitude for 7. According to the 
Landau-de Gennes theory, the rescaled bulk nematic order parameter at the phase transition 
is Sb = 0.61 which corresponds to the measured typical value of S m b — 0.3 for the Maier- 



6 



EUROPHYSICS LETTERS 




Fig. 3 - a) Three-dimensional phase diagram of the different wetting layers for 300 nm particles. 
Surface with solid grid lines: transition uniaxial oblate-boojum, surface with dashed grid lines (hardly 
visible): first-order transition boojum-tetrahedral. The solid line on the upper surface is a tricritical 
line. It separates a first-order (smaller 7) from a second-order (larger 7) transition. Note the very 
narrow existence region of the boojum. b) Contour plots of the three-dimensional phase diagram for 
various So. 



Saupe order parameter S m — (3 cos 8 — l)/2 [20]. For a perfectly aligned nematic (prolate 
order) S m = 1; the perfect oblate order gives S m = —1/2. This means that in our scaling 
the minimum value of So in the surface potential is So = — Sb/S m b/2 = —1.02. This value 
is further decreased for smaller S m b, so that our predicted transitions from the oblate to the 
biaxial surface phases, especially the tetrahedral configuration, should be observable. 

Finally, we add some comments about the influence of fluctuations on our phase diagrams. 
It was argued that in the planar geometry, the transition from the biaxial to the oblate 
wetting layer is governed by defect unbinding, i.e., by the Berezinskii-Kosterlitz-Thouless 
(BKT) mechanism [11-13]. As a result, the second-order transition line is shifted to smaller 
temperatures (see e.g. Ref. [12]). The treatment of fluctuations is beyond the scope of this 
paper. However, it seems possible that under the BKT mechanism the boojum configuration 
vanishes completely from the phase diagrams in Fig. and H due to its very narrow existence 
region. In an elaborate and elegant treatment, Nelson finds that the tetrahedral configuration 
of the two-dimensional nematic is stable against thermal fluctuations deep in the nematic 
phase, whereas close to the phase transition the fluctuations are fairly large [1]. How these 
results apply to our studies is not completely clear. It certainly depends on how strong the 
biaxial ordering is developed in the tetrahedral configuration. 




Fig. 4 - The same phase diagrams as in Fig. but now for 100 nm particles. 
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In conclusion, we have shown that nematic wetting layers around spherical particles possess 
indeed a tetrahedral configuration where 1/2 disclinations are located at the vertices of a 
tetrahedron. This offers an interesting route towards tetravalent colloids. Planar anchoring, 
necessary for the realization of this tetrahedral configuration, can be realized in inverted 
nematic emulsions [9] , or in lyotropic liquid crystals conmposed of rodlike surfactant micelles 
[21]. A possible method to detect nonuniform wetting layers is depolarized light scattering 
[22]. Furthermore, the Brownian rotational motion of the particles should be measurable by 
dynamic light scattering [23]. 

* * * 
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